An astronomical search for evidence of new physics: Limits on gravity-induced 
birefringence from the magnetic white dwarf RE J031 7-853 
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The coupling of the electromagnetic field directly with gravitational gauge fields leads to new 
physical effects that can be tested using astronomical data. Here we consider a particular case for 
closer scrutiny, a specific nonminimal coupling of torsion to electromagnetism, which enters into 
a metric-affine geometry of space-time. We show that under the assumption of this nonminimal 
coupling, spacetime is birefringent in the presence of such a gravitational field. This leads to the 
depolarization of light emitted from extended astrophysical sources. We use polarimetric data of 
the magnetic white dwarf RE J03f7-853 to set for the very first time constraints on the essential 
coupling constant for this effect, giving k 2 < (19 m) 2 . 

PACS numbers: 04.80.Cc 
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I. INTRODUCTION 

The designation new physics is used here to refer to 
phenomena beyond the scope of the standard model of 
particle physics or of general relativity. Today, efforts to 
develop a quantum theory of gravity or a complete, con- 
sistent and unified theory of matter and all its interac- 
tions are rich and compelling sources of speculation about 
new physics. Certainly, the effective field theories that 
emerge as low-energy limits of string theories are littered 
with new fields and interactions. Since cherished sym- 
metries like CPT and Lorentz invariance can be broken 
in these contexts 0, , high-precision experimental and 
observational tests of these symmetries offer particularly 
effective ways of searching for evidence of new physics. 

In this letter we focus on consequences of interactions 
that couple the electromagnetic field to new physics. In 
recent papers by Kostelecky & Mewes Q, couplings of 
the electromagnetic field to an axial-vector, {Uaf) 11 -, an d 
to a tensor, (kp) fJ -'" 7T , are used to provide a generic char- 
acterization of conceivable new physics influencing the 
photon sector. The Uaf and kp structures could, for 
example, be constructs of vector and tensor fields in an 
underlying theory. If such fields acquire nonzero vacuum 
expectation values, their couplings to the electromagnetic 
field spontaneously break global Lorentz invariance. The 
kAF structure breaks CPT as well. While motivated 
by current speculations regarding fundamental physics, 
the conceptual framework and formalism of the work by 
Kostelecky and Mewes is that of Lagrangian-based local 
field theory. With different motivations, questions re- 
garding the most general way in which new fields may 
couple to the electromagnetic field have been addressed 
in this framework before. It is reassuring to find all of 
these studies reaching equivalent conclusions. 

Today's strong empirical foundation for general rel- 
ativity is a consequence of continuing efforts to probe 



for evidence of this sort of new or alternative physics, 
see Will Q for a review. In this research setting the 
Lagrangian-based local field theory framework is often 
called the Dicke framework Q. Theoretical alternatives 
to general relativity can feature gravitational fields in ad- 
dition to the symmetric second-rank tensor potential of 
general relativity. One distinguishes two classes of alter- 
native theories based on the way in which these gravi- 
tational fields couple to matter. A theory is metric if it 
admits a representation in which matter fields are all min- 
imally coupled in the familiar way to a symmetric second- 
rank tensor gravitational field alone, implying that aux- 
iliary fields couple only to that metric tensor field and to 
each other. Otherwise, a theory is nonmetric. It is in the 
study of nonmetric theories which couple gravitational 
fields directly to the electromagnetic field that one finds 
analogues of the recent work by Kostelecky and Mewes. 
The equivalent of their Lagrangian density governing the 
dynamics of the electromagnetic field can be found in 
Ni's X9 formalism p|. There, this Lagrangian density is 
expressed in terms of a fourth-rank tensor field, x"^** 6 > 
constructed from a symmetric second-rank tensor field, 
g^ u , and auxiliary gravitational fields, 

Cem = X a0l5 F a0 F jS . (1) 

Ni noted that theories encompassed by this formalism 
can predict birefringence and used pulsar polarization ob- 
servations to constrain this possibility twenty years ago. 

Recognition that gravitation theories touted as viable 
alternatives to general relativity could predict such bire- 
fringence stimulated new interest in such possibilities in 
the 1990s and sharp constraints were imposed on birefrin- 
gence predicted by Moffat's non-symmetric gravitation 
theory (NGT) 0, II H E3, EH ■ The analysis of birefrin- 
gence by Kostelecky and Mewes and that performed 
by Haugan and Kauffmann ^lj m the general context 
provided by the \9 formalism run parallel to one another. 
Indeed, the 3+1 decomposition of the x tensor made by 
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Haugan and Kauffrnann yields SO (3) tensors equivalent 
to the k matrices of Kostelecky and Mewes, 

{K DE ) ij = 2C j , (khbY j = 2C J , {kdbY 3 = H 3 - (2) 

Kostelecky and Mewes restrict their attention to cases 
in which these quantities arc spacetime constants. Hau- 
gan and Kauffrnann allow them to vary, reflecting the 
structure of fields involved in new physics that are gener- 
ated by localized sources. Their work establishes that 
evidence of such birefringence can reveal new physics 
even when the fields involved do not acquire nonzero vac- 
uum expectation values and global Lorentz invariance is 
not spontaneously broken. String-inspired speculation 
regarding new physics provides a compelling motivation 
to continue searching for evidence of birefringence of the 
kind analyzed in preceding references. 

As additional motivation we offer a new example, 
drawn from the gravitation physics literature, of a the- 
oretical alternative to general relativity that naturally 
predicts such birefringence if we assume a special non- 
minimal coupling of torsion to electromagnetism. Be- 
sides Moffat's NGT this provides a new concrete test case 
against which one can gauge the precision of birefringence 
searches. In this context it is important to note that the 
new version of Moffat's NGT |l2j also provides a valid 
alternative to general relativity which has overcome the 
technical difficulties found by Damour et al. and 
Clayton et al. 0]. The foil for general relativity we 
present in section II, is drawn from the class of metric- 
affine theories of gravity . We hope that this exam- 
ple will prompt further work on gravity-induced birefrin- 
gence predicted by metric-affine theories and other alter- 
natives to general relativity; see In this paper we 
also report new observational constraints on such gravity- 
induced birefringence based on data from white dwarfs - 
which constitute the most abundant stage of final stellar 
evolution - possessing a strong magnetic field. 



II. GRAVITY-INDUCED BIREFRINGENCE IN 
METRIC-AFFINE THEORIES 

Members of the class of metric-affine theories of gravity 
feature torsion and/or nonmetricity gravitational fields in 
addition to a symmetric second-rank tensor gravitational 
potential |l5j . The minimal or geometrical coupling 
scheme for coupling gravity to matter in metric-affine 
theories does couple torsion and nonmetricity to some 
matter fields so metric-affine theories are, in fact, non- 
metric theories. However, the minimal coupling scheme 
does not couple torsion or nonmetricity to the electro- 
magnetic field |18| . Only nonminimally coupled metric- 
affine theories can conceivably predict birefringence. 

Here, we focus entirely on the effect of nonminimal cou- 
pling to torsion. We suggest, in addition to the conven- 
tional Maxwell Lagrangian, the additional nonminimal 
piece |2(J 



see 0] for other possibilities, where k is a coupling con- 
stant with the dimension of length, * denotes the Hodge 
dual, T denotes the torsion and F the electromagnetic 
field, which is related in the usual way to its potential A. 
This addition is consistent with gauge invariance and, so, 
with charge conservation. It corresponds to a term 

SCem =Sx al37S F a pF 7S , (4) 

in the notation of Haugan and Kauffrnann 11]. This 
correspondence permits us to infer consequences of the 
nonminimal torsion coupling above from the results of 
the general analysis of reference |Tlj . Since our immedi- 
ate interest is the interpretation of stellar data, we focus 
on the case of static, spherically symmetric metric-affine 
fields to explore the dominant effects of solar torsion on 
the propagation of light from the Sun. Tresguerres 0] 
has expressed the most general form of a static, spheri- 
cally symmetric torsion field in terms of six radial func- 
tions. In his notation these are a, f3 and 7(j) for i =1 to 
4. Plugging his representation into SCem above yields 
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where E and B refer in the usual way to the electric and 
ma gnet ic components of F in Tresguerres' coordinates 
|26l l27j . Note that the three- vector index 1 refers to the 
radial direction. In terms of the components of the £, 
£ and 7 tensors of reference |Tlj . we can read off from 
the corresponding expression for SCem bow the nonzero 
components of the tensors £, £ and 7 are expressed in 
terms of the Tresguerres functions for the nonminimal 
coupling considered here. To make a specific prediction 
of the effect this torsion coupling will have on the prop- 
agation of light from a star, we must consider a specific 
metric-affine candidate for the stellar field. We choose a 
special case of a static, spherically symmetric field found 
by Tresguerres |20| in which nonmetricity vanishes, cor- 
responding to vanishing dilatation and shear charge. For 
simplicity we also set the cosmological constant to zero. 
The corresponding forms of Tresguerres' a, (3 and 7m 
functions are 
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1/1 m\ , . 1/1 m. 

-9 - + I2 >7(3)(r) = r --3 (7) 



with 7(2) = 7(4) = 0. A somewhat tedious recapitulation 
of the analysis in ^lj in this case leads to the conclusion 
that, to leading order, the specified nonminimal torsion 
coupling and chosen stellar torsion field singles out linear 
polarizations with a fractional difference in their propa- 
gation speeds of 
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where 9 denotes the angle relative the outward radial di- 
rection in which the light is propagating. In this context 
it is interesting to compare (JHJ with equation (21) and 
(22) of Eft 

We are interested in the effect that this differential 
propagation has on the polarization of light as it trav- 
els from a localized source on the stellar surface to an 
observer. This is determined by the phase shift that ac- 
cumulates between the polarization components singled 
out by the stellar field as the light propagates. Recapit- 
ulating the analysis performed in [3, llOfl , we find that in 
this case 
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where ^ denotes the cosine of the angle between the line 
of sight and the normal on the stellar surface (/i = 1: stel- 
lar disk center, [i = 0: limb), A is the light's wavelength 
and R is the stellar radius. 



III. DATA, ANALYSIS AND RESULTS 

All light received from a pointlike source of polarized 
radiation in a gravitationally birefringent environment 
suffers the same phase shift A$(/x). Introducing wave- 
length dependent Stokes parameters I\, Q\, U\,V\ to 
describe polarized light [21|. with Stokes Q defined to 
represent the difference between linear polarization par- 
allel and perpendicular to the local stellar limb, gravi- 
tational birefringence introduces a crosstalk between the 
linear polarization parameter Stokes U and the net circu- 
lar polarization, V. This crosstalk is such that although 
the observed values U i, s and V^s differ from the values 
emitted by the source, U SIC and V STC , the composite de- 
gree of polarization remains equal: {U 2 hs + V 2 ^) 1 ^ 2 = 
([/ s 2 rc + K, 2 .,,) 1 / 2 . If an extended source covering a range 
of (i values is observed then light emitted from different 
points suffers different phase shifts and, so, adds up to 
an incoherent superposition. Summing over the differ- 
ent contributions, while using the additive properties of 
Stokes parameters, yields a reduction of the observed po- 
larization relative to the light emitted from the source: 



(^bs 



) 1/2 < (U 2 TC + U s 2 rc ) 1/2 . Since the rotation- 



ally modulated polarization from magnetic white dwarfs 
can only be produced by an extended source any 
observed (i.e. non-zero) degree of polarization provides 
a limit on the strength of birefringence induced by the 
star's gravitational field @. 

The polarized radiation from white dwarfs is produced 
at the stellar surface as a result of the presence of ultra- 
strong (up to 10 5 T) magnetic fields [23. Since the disk 
of a white dwarf is unresolved, only the total polarization 
from all surface elements is observable: 



VxMk 2 ) = ^J J V x (ti,B,B\\) cos(A$) M c$# 



(10) 



where the Stokes parameter V\ changes over the visible 
hemisphere and depends on the wavelength A, the loca- 
tion /j, (limb darkening), the total magnetic field strength 
B(8, (f>), the line-of-sight component B\\ (8, <fr), and on the 
parameters of the stellar atmosphere. Gravitational bire- 
fringence reduces the polarization by means of cos(A<i>). 
Theoretically, the Stokes parameters can be calculated by 
solving the radiative transfer equations through a magne- 
tized stellar atmosphere on a large number of surface ele- 
ments on the visible hemisphere (e.g. H^)- If the star is 
rotating, the spectrum and polarization pattern changes 
according to the respective magnetic field distribution 
visible at a particular moment. To obtain the degree of 
circular polarization, Eq. (|l(Jfl has to be divided by the 
total stellar flux 2a. tot emitted to the observer at wave- 
length A. Below we will calculate a maximum circular 
polarization VA. max //A.tot from radiative transfer calcu- 
lations which is in general higher than the observed value 
^A,obs/-^A,obs- Then we assume that the reduction from 
VA,max to V\ i0 bs i s entirely due to the factor cos(A<I>(£; 2 )) 
in Eq. l|lUfl thereby calculating the maximum value for 
k 2 , i.e. our limit on k 2 is reached as soon as V^tot/^tot 
in Eq. I|10[l becomes smaller than I^A.obs/A.tot for a cer- 
tain value of k 2 . 

RE J0317-853 is a remarkable object within the class of 
isolated magnetic white dwarfs. Besides being the most 
rapidly rotating star (P = 725 sec) of this type it is also 
the most massive at 1.35 M©, close to the Chandrasekhar 
limit [2^|. The corresponding radius is only 0.0035 2?©. 
With a reported VA i0 bs/2A,tot °f 20% at A = 576 nm 
31], RE J0317-853 is also the magnetic white dwarf with 
the highest known level of circular polarization. Due to 
its small radius and high degree of circular polarization, 
RE J0317-853 is a very suitable object for setting limits 
on gravitational birefringence. The analysis of time re- 
solved UV flux spectra obtained with the Hubble Space 
telescope has shown that the distribution of the field 
moduli is approximately that of an off-centered magnetic 
dipole oriented obliquely to the rotation axis with a polar 
field strength of Bd = 3.63 • 10 4 T, leading to visible sur- 
face field strengths between 1.4-10 4 T and 7.3-10 4 T [3(|. 
This model is not only able to describe the UV, but also 
the optical spectra (Jordan et al. in prep.), which means 
that the distribution of the magnetic field moduli - but 
not necessaryly of the longitudinal components - is cor- 
rectly described. This result is completely independent 
of the magnitude of the gravitational birefringence. From 
radiative transfer calculations it follows that at the phase 
of rotation when the maximum value of 20% polarization 
at 576 nm is measured, almost the entire visible stellar 
surface is covered by magnetic fields between 1.4- 10 4 and 
2.0 • 10 T, with only a small tail extending to maximum 
field strengths of 5.3 • 10 4 . This distribution can best be 
modeled by assuming a rotational phase where the axis of 
the off-centered dipole is nearly perpendicular to the line 
of sight. Using this special field geometry we calculated 
a histogram distribution of the visible surface magnetic 
field strengths in order to set sharp limits on gravita- 
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tional birefringence. For each field strength bin of the 
histogram we calculated the maximum circular polariza- 
tion from radiative transfer calculations by assuming that 
the field vector always points towards the observer. The 
total maximum polarization from the whole visible stel- 
lar disk without gravitational birefringence is then cal- 
culated by adding up the contributions from each field 
strength bin weighted with its relative frequency. This 
results in Vx,max/ I\,tot = 26.5%. Assuming that the re- 
duction to VA,obs/^A,tot = 20% is entirely due to gravity 
induced depolarization - and not due to the fact that in 
reality not all field vectors point towards the observer - 
we find an upper limit for this effect of k 2 < (19 m) 2 . Since 
there is always a small uncertainty in determing the ex- 
act mass of a white dwarf, we also calculated an upper 
limit on k 2 assuming a lower mass of 1 M Q . This leads to 
k 2 <(22m) 2 . An even more extreme assumption would 
be to use the maximum circular polarization predicted 
by all field strengths in the interval 1.4 • 10 — 5.3 ■ 10 4 
(reached at 5.3 • 10 4 T) for the whole stellar disk. Then 
we obtain VA,max/-T\,tot = 48.3% and an upper limit of 
k 2 <(30.5 m) 2 . Independent from any dipole model and 
without any reference to radiative transfer calculations 
the assumption of 100% emerging polarization leads to 
fc 2 <(45m) 2 . In order to compare the quality of our 
method with previous similar results from Solanki et al. 
Q we also set new upper limits on the NGT parame- 
ter £ 2 wich also causes gravitational birefringence in case 
of £ 2 7^ 0. Assuming VA ima x/A,tot = 26.5% leads to 
£ 2 <(1.8km) 2 in contrast to ^ 2 <(4.9km) 2 , determined 
for the white dwarf GRW +70°8247. 



IV. CONCLUSIONS 



We have used spectropolarimetric observations of the 
massive (1.35 Mq) white dwarf RE J0317-853, to impose 
new strong constraints on the birefringence of space- 
time in the presence of a gravitational field. Such an 
effect is predicted by new physics in the form of a direct 
coupling of the electromagnetic field to auxiliary grav- 
itational fields. Since the gravity-induced birefringence 
of space violates the Einstein equivalence principle, our 
analysis also provides a test of this foundation of gen- 
eral relativity and other metric theories of gravity. We 
consider as a specific case a metric- affine theory that cou- 
ples the electromagnetic field with torsion and for which 
a static spherically symmetric field has been found |2Cj . 
The data provide an upper limit for the relevant cou- 
pling constant k 2 of (19 m) 2 - or (45 m) 2 for the most 
conservative assumptions on. Considerably tighter limits 
based on the same astronomical source could be provided 
by measurements of circular polarization in the FUV (in 
particular associated with Lya absorption features) and 
also by a consistent model for the magnetic field geometry 
which reproduces the spectropolarimetry measurements 
in the optical. 

We are grateful to Dayal T. Wickramasinghe and 
Friedrich W. Hehl for helpful remarks and valuable dis- 
cussions. This research has made use of NASA 's Astro- 
physics Data System Abstract Service. 
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